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ABSTRACT 


The problem of electromagnetic (EM) scattering from irregularly shaped, thin, metallic flat 
plates in free space is solved using the electric field integral equation (EFIE) approach in 
conjunction with the method of moments (MoM) with quadrilateral meshing. An irregularly 
shaped thin plate is discretized into quadrilateral patches and the unknown electric surface 
current over the plate is expressed in terms of proper basis functions over these patches. The 
basis functions for the electric surface current density that satisfy the proper boundary conditions 
on these quadrilateral patches are derived. The unknown surface current density on these 
quadrilateral patches is determined by setting up and solving the electric field integral equation 
by the application of the MoM. From the knowledge of the surface current density, the EM 
scattering from various irregularly shaped plates is determined and compared with the earlier 
published results. The novelty in the present approach is the use of quadrilateral patches instead 
of well known and often used triangular patches. The numerical results obtained using the 
quadrilateral patches compare favorably with measured results. 

1. INTRODUCTION 

The EM scattering from a complex-shaped metallic object is of practical interest to 
electromagnetic analysts as well as to engineers. The EM scattering from polygonal shaped 
plates is of special interest because a complex-shaped object can always be modeled as an 
interconnection of these polygonal flat plates. An EM scattering analysis of polygonal shaped 
metallic plates can be accomplished by using various numerical techniques [1-4]. One of the 
widely used techniques is the MoM in which a polygonal plate is first discretized into number of 
triangular patches called sub domains. The unknown surface current density on these sub 
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domains is then expressed in terms of well known Rao-Wilton-Glisson (RWG) [3] basis 
functions. The unknown surface current density on these sub domains is then determined using 
the EFIE in conjunction with the MoM. Although the triangular meshing and use of RWG basis 
functions is a very popular and mature concept, the non-vanishing normal component of surface 
current on one of the inclined edges of sub-domain triangles attached to the open /boundary edge 
of polygonal plate will give a non-zero normal component of the surface current. However, the 
non-zero normal component of the surface current on the boundary edge will approach zero 
when very fine discretization is used in the region close to the open/boundary edges. This 
problem will not arise if one uses a sub domain element (near the boundary) whose edges will be 
either along the boundary or normal to the boundary (such as quadrilateral element). Another 
approach as implemented in [4] divides a polygonal plate into small rectangles. The surface 
current density over each sub domain is expressed in terms of overlapping triangular functions 
(roof-top functions) in the direction of current flow and a pulse function in the orthogonal 
direction. The EFIE in conjunction with the MoM is then used to determine the surface current 
density. An important advantage of this procedure is that the impedance matrix encountered in 
this procedure has a block Toeplitz nature, a property that is useful in reducing the matrix filling 
time. However, the major drawback of this scheme is that the curve boundaries are 
approximated by stair cases and hence for accurate results requires very fine discretization. 

In this report the problem of electromagnetic scattering of plane waves by an arbitrarily 
shaped, thin metallic plate is studied using the EFIE method in conjunction with quadrilateral 
meshing. Use of the equivalence principle and the free space Green’s function in deriving the 
EFIE with the surface current density as an unknown is described in detail. With proper choice 
of the expansion and testing functions over the quadrilateral mesh, the MoM is described to 
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convert the EFIE into a matrix equation. The procedure to obtain the scattered far field by the 
plate in terms of the electric surface current is also described. 

The remainder of report is organized as follows. The formulation of the problem in 
terms of the EFIE using the surface equivalence principle is developed in section 2. Numerical 
results on the mono-static radar cross section of various polygonal shaped metallic plates are 
presented in section 3 and compared with earlier published data. The advantages and limitations 
of the present formulation are discussed in section 4. 

2. THEORY 

Consider a time harmonic electromagnetic plane wave incident on an irregularly shaped, 
infinitesimally thin metallic plate as shown in Figure 1 . The incident field with a time 
dependence may be written as 

Incident 



Figure 1 : Geometry of arbitrarily shaped thin metallic plate in x-y plane 
and illuminated by a plane wave 
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( 1 ) 


E„^{dE,,+(t>^E,^)e 




where = -k^ cos(<z>, ) sin(^,. ) + y sin(<z),. ) sin(^, ) + £ cos(^, )] , r = + yy + & , 

Eg =|£,|cos(oro), and =|e, sin(ao), being the free- space wave number. With 

reference to Figure 1, cto 0 corresponds to H-polarization and = 90® corresponds to E- 
polarization. From equation (1), the x-, y-, and z-components of the incident field may be 
written, respectively, as 

E^. = cos(6>, ) cos(<z); ) - sin(^z), ) (2) 


£v, = cos(^,.)sin(<z),.)-£^ cos(^,.) 

£.. =-Eg^ sin(6>, ) 

The corresponding magnetic field components are obtained through 

1 r 


(3) 

(4) 


W.., = 






(5) 


where //g is the free-space impedance. The incident field with £^ Oand £^ = 0 is called the 
H-polarized wave and £^ = 0 and £^ 0 is called the E-polarized wave. 

To facilitate the solution of the problem using the equivalence principle, the thin 
metallic plate is replaced by the equivalent electric surface current density J ( jt, y) . The 
scattered field due to J may be obtained through the vector magnetic potentials as 

E^{J) = -jaA-V^ (6) 
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( 7 ) 


H^(J) = —VxA 
Mo 

where the magnetic vector and scalar potentials are given by 


-Jkolr-i"] 


4;r ' r - r 




^=^r|(v.jr)£^ 

4;r<y£’o ^ \r-r 


-ds 


( 8 ) 

(9) 


where £q and jUq are the permitivity and permeability, respectively, of free-space, f and r' are 

respectively the coordinates of field and source points. From the knowledge of the scattered 
field and incident field, the EFIE is set up by equating the total tangential electric field over the 
surface of the polygonal plate to zero. Hence 


b<!rt+v«iL=[£.L (10) 

To determine the approximate solution of the integral equation in (10) using the MoM, 
the polygonal surface is approximated by a union of quadrilaterals as shown in Figure 2. On the 
n"" quadrilateral the electric surface current density can be represented by a superposition of 
vector functions B,„. (a:, y) as 

J,(.x.y) = j^Tj„(x.y) ( 11 ) 

/=1 

where is the amplitude of electric current normal to the i''' edge of the n"' element 
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Node #1 


Figure 2: Arbitrarily shaped thin metallic plate discretized into 
quadrilaterals, and are typical elements with their four nodes 


(as shown in Figure 4). The vector basis function (jc, y) representing the surface current 
density on the n"' quadrilateral is required to satisfy the appropriate boundary conditions; such 
as the normal component of y) to the i"’ edge must be unity, and must have a value of 

zero on the other remaining edges. Construction of such a vector basis function with help of 
Figure 3 is explained in detail in the Appendix A. 


(X3. ys) 
Node #3 



T| 

A 


No(de #2 



Node #1 


Edge #4 

n=-i 


Node #4 


Figure 3; (a) n''’ quadrilateral in x-y plane, (b) n^'' quadrilateral in if plane 
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Node #6 



Figure 4: (a) quadrilateral with four edges, (b) n^^and its adjacent quadrilateral 





( 12 ) 


where A,„. and are obtained from equations (8) and (9), respectively, by replacing J by 


5,„- . Selecting B^j with j = 1,2, 3, 4 and m = 1,2,3,.. .N as a testing function and application of 
the MoM to equation (12) yields 




IZr. 

1=1 


Ajt 


■JJb.A.JJb 


■ds^ + 






AjtCOEr 


J J jv • ds„, = - j J • [£,„ J ds„ 

J m tan 


1 


(13) 


The left side of equation (13) over the n'* quadrilateral (Figure 4(a) ) can be written as 



Ati.n 
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( 14 ) 
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■ds„ + 




A;tctje, 


n m 'I 


( 15 ) 


Now consider union of two adjacent quadrilaterals with a common edge as shown in Figure 
4(b). The left side of equation (13) for two adjacent quadrilaterals can be written as 
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In the above expression, the first letter of superscript represents the test quadrilateral and the 
second letter represents the source quadrilateral element. To ensure continuity of current at the 
common edge it is assumed that ^ . With an assembled mesh of all quadrilaterals 

over the surface of polygonal plate, the left side of equation (13) can be written as 

[z][r] (16) 
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where the elements of the matrix are obtained from A""' and the column matrix [r ] represents 

the unknown amplitude of current at the edges of quadrilateral elements. Using (16) the equation 
(13) can be written in a matrix form 

[z][r] = [y] (17) 

where the elements of column vector matrix are obtained from 

ds. (18) 

m tan 


The edge currents obtained after solving the matrix equation (17) are then used to determine the 
EM scattering from polygonal thin metallic plates. It is noted that the 7^. on the boundary edges 
are forced to zero before solving the matrix equation (17). 

Scattered Far Field: 

The scattered far field due to the polygonal plate is estimated from 


^6 = -jKVoiK cos(<z>)cos(^) + sin(</))cos(<9)) (19) 

= -JKno (“Av sin(<Z>) + A^ cos(<p)) (20) 

where A^ and A^ are obtained from 


,-jkor N 4 


A = - 


4/r-ri3'1- 


^ jkQ ( a' cos( ^ ) sin( ^ )+ y ' sin( ^ ) sin( ^ 


( 21 ) 


From the far fields, the radar cross sections, for co and cross polarizations are obtained as 
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( 22 ) 


(23) 

(24) 

(25) 

3. NUMERICAL RESULTS 

In this section the Radar Cross Section (RCS) of polygonal, thin metallic plates of various 
shapes is calculated using the quadrilateral meshing scheme. To validate the present analysis 
and the computer code developed, the numerical results obtained by the present method are 
compared with (1) the RCS computed using the ESP code [4], (2) the RCS computed using 
triangular meshing', and (3) the RCS measured in the Langley Experimental Test Range Facility 
[4]. 

RCS of Hexagonal Plate: 

As a first example, a hexagonal plate with a = 2.074 cm, as shown in Figure 5, is 
considered. Using the hexa.SES file (listing given on the attached disk), the hexagonal plate is 
modeled and meshed (also shown in Figure 5 using the Geostar/COSMOS commercial 
software). The hexa.MOD file generated by Geostar is then is imported into PoMePl (Polygonal 
Metallic Plate, listing given on the attached disk ) electromagnetic code written in FORTRAN to 
estimate the RCS of polygonal, thin metallic plates. The mono-static RCS of the hexagonal plate 

' MoM code developed at the Electromagnetic Research Branch of NASA Langley Research Center by the author. 


^hh 
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^EH ^ 


^<1 


7— 77 

\E 
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is calculated for both E- and H-polarized incident waves for various discretization levels and is 
shown in Figures 6-7. It is clear from the Figures 6-7 that the discretization level with the 
number of edges more than 232 yields fairly consistent numerical results. Hence for the RCS 
computation of hexagonal plate shown in Figures 8-11 the number of edges used was 232. The 
Figures 8-1 1 also give a comparison of the results obtained using the quadrilateral meshing with 
the results obtained using triangular meshing, the ESP code [4], and the measured data [4 ]. 
From Figures 8-1 1 it can be observed that the results obtained using the quadrilateral meshing 
agrees well with the other numerical data as well as with the measured results. However, for E- 



Figure 5: Thin, metallic flat plate of hexagonal shape with side a = 2.074 
cm and descritized into quadrilateral sub-domains. 
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Figure 6: Monostatic RCS of hexagonal plate shown in Figure 5 when 
illuminated by E-polarized plane wave as a function of number of edges 



Figure 7: Monostatic RCS of hexagonal plate shown in Figure 5 when 
illuminated by H-polarized plane wave as a function of number of edges 
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Figure 8 Monostatic RCS of hexagonal plate (shown in Figure 5) when illuminated by E-polarized plane wave 
OCq = 90*’ at = —90" and = 90*^, solid line- present method, dashed line RWG functions, solid square 
symbol- measurement, frequency = 1 1.81 1 GHz. 
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Figure 9 Monostatic RCS of hexagonal plate (shown in Figure 5) when illuminated by H-polarized plane 
wave dfo = 0" at = —90" and (f>j = 90“ , solid line- present method, dashed dot line ESP-code, solid 
square symbol- measurement, and dashed line RWG function, frequency = 1 1.81 1 GHz. 


13 



Figiire 10 Monostatic RCS of hexagonal plate (shown in Figure 5) when illuminated by H-polarized 
plane wave = 0" at = 0" and (p^ = 180®, solid line- present method, dashed dot line ESP- 
code, solid square symbol- measurement, and das ed line RWG function, frequency = 1 1.81 IGHz. 

-m~ 



^,,deg. ► 

Figure 1 1 Monostatic RCS of hexagonal plate (shown in Figure 5) when illuminated by E-polarized 
plane wave (Xq = 90" at = 0" and (p. = 180®, solid line- present method, dashed line RWG 
functions, solid square symbol- measurement, frequency = 1 1.81 IGHz 
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polarized incidence close to grazing angles, the computed results do not agree well with the 
measured data. This is due to the fact that in the computation the finite thickness of the plate is 
assumed to be zero. 

RCS of Equilateral Triangular Plate: 

A equilateral triangular plate as shown in Figure 12 is considered as a second example 
for the validation of the present code based on the quadrilateral meshing. Using the 
COSMOS/GEOSTAR the equilateral triangle is discretized into quadrilateral patches as shown 
in Figure 12. The *.MOD file, in conjunction with PoMePl code is used to estimate the 
monostatic RCS of the equilateral triangular plate. The results of the computations are shown in 
Figures 13-18 along with the experimental and other numerical data. The monostatic RCS 



Figure 12 Equilateral triangular, thin, metallic flat plate with a = 5.08 cm 

calculated using the quadrilateral meshing compares well with the results obtained using the 
triangular meshing and the ESP - code. The experimental data also compares well with the 
numerical data as long as the angle of incidence is not close to grazing. The discrepancy 
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between the measured and numerical data at angles close to grazing is because of the 
assumption of zero thickness in the numerical simulation. The discrepancy between measured 
and numerical data as a function of frequency for the grazing incident angles, as shown in Figure 
17, illustrate the necessity of taking into account the thickness of the plate in the numerical 
simulation techniques for more accurate results. 
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Figure 13 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by E-polarized plane wave = 90" at <z>, = 0" and =180°, solid line- present 
method, short dashed line- RWG functions, long dashed line ESP-code, solid square symbol- 
measurement, frequency = 1 1.81 IGHz. 
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Figure 14 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by H-polarized plane wave = 0" at = 0° and = 180°, solid line- present 

method, short dashed line- RWG functions, dashed dot line ESP-code, solid square symbol- 
measurement, frequency = 1 1.81 1 GHz 



Figure 15 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by E-polarized plane wave = 90" at = -90" and = 90° , solid line- 


present method, short dashed line- RWG functions, dashed dot line ESP-code, frequency = 
11.811GHz 
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Figure 16 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by H-polarized plane wave = 0" at = -90" and </>i = 90*^ , solid line- 
present method, short dashed line- RWG functions, dashed dot line ESP-code, solid square 
symbol- measurement , frequency = 1 1.81 1 GHz. 
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Figure 17 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by E-polarized plane wave = 90" as a function of for frequency = 

1 1.81 1 GHz, 6. = 90" . solid line- present method, short dashed line- RWG functions, 
dashed dot line ESP-code, frequency = 1 1 .81 IGHz. 
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Figure 18 Monostatic RCS of equilateral triangular plate (shown in Figure 12) when 
illuminated by E-polarized plane wave = 90'’ as a function of frequency for 

(!>j = 90 and = 90" . solid line- present method, short dashed line- RWG functions, 
dashed dot line ESP-code, solid square symbols- measurement. 


RCS of Equilateral Triangular Plate With Concentric, Equilateral Triangular Hole. 

An equilateral triangular plate with a concentric equilateral triangular hole as shown in 
Figure 18 is considered as an another example for the validation of the present code based on the 
quadrilateral meshing. Using the COSMOS/GEOSTAR the equilateral triangle is discretized 
into quadrilateral patches as shown in Figure 18. The *.MOD file, in conjunction with PoMePl 
code is used to estimate the mono-static RCS of the equilateral triangular plate with an 
equilateral triangular hole. The results of the computations are shown in Figures 19-24 along 
with the experimental and other numerical data. The mono-static RCS calculated using the 
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quadrilateral meshing compares well with the other numerical data as well as with the 
experimental results. 



Figure 19 Equilateral triangular thin metallic flat plate with a concentric equilateral 
triangular hole, a = 5.08 cm, b = 2.54 cm and descritized into quadrilateral patches. 

For this example, it is also observed that the numerical simulation data do not match well 
with the experimental results for grazing incidence. 


20 






Figure 20 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown in 
Figure 19) when illuminated by E-polarized plane wave CCq = 90*’ at = —90° and = 90^ , solid 

line- present method, short dashed line- RWG functions, dashed dot line ESP-code, solid square 
symbols- measurement, frequency = 11.811GHz. 



Figure 21 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown in 
Figure 19) when illuminated by H-polarized plane wave CX^ = 0° at = —90° and = 90** , solid 

line- present method, short dashed line- RWG functions, dashed dot line ESP-code, solid square 
symbols- measurement, frequency = 11.811GHz. 
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Figure 22 Monostatic RCS of equilateral triangular plate with equilateral triangular hole 
(shown in Figure 19) when illuminated by E-polarized plane wave OTq = 90" at (p. = 0" and 

(pj = 180** , solid line- present method, short dashed line- RWG functions, dashed dot line ESP- 
code, solid square symbols- measurement, frequency = 1 1.81 IGHz. 
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Figure 23 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown 
in Figure 19) when illuminated by H-polarized plane wave - 0" at (Zi, = 0" and (p^ = 180**, 

solid line- present method, short dashed line- RWG functions, dashed dot line ESP-code, solid 
square symbols- measurement, frequency = 11.811GHz. 
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Figure 24 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown in Figure 19) 
when illuminated by E-polarized plane wave cx^ =90" as a function of for frequency = 1 1.81 1 GHz, 

6j = 90" . solid line- present method, short dashed line- RWG functions, dashed dot line ESP-code, solid 
square symbols- measurement, frequency = 11.811GHz. 
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Figure 25 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown in Figure 19) 
when illuminated by E-polarized plane wave =90" as a function of frequency for (Zi, = -90 and 

6j = 90" (grazing incidence at edge), solid line- present method, short dashed line- RWG functions, dashed 
dot line ESP-code, solid square symbols- measurement 
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Figure 26 Monostatic RCS of equilateral triangular plate with equilateral triangular hole (shown 
in Figure 19) when illuminated by E-polarized plane wave = 90" as a function of frequency 

for = 90 and 0 ^ — 90" (grazing incidence at tip), solid line- present method, short dashed line- 
RWG functions, dashed dot line ESP-code, solid square symbols- measurement 

RCS of Diamond Shaped Plate 

A diamond shaped polygonal plate as shown in Figure 25 is considered as a last example 
for validation of the present code. Using the COSMOS/GEOSTAR the diamond shaped plate is 
discretized into quadrilateral patches as shown in Figure 25. The *.MOD file, in conjunction 
with PoMePl code is used to estimate the mono-static RCS of the diamond shaped plate. The 
results of the computations are shown in Figures 26-30 along with the experimental and other 
numerical data. The mono-static RCS calculated using the quadrilateral meshing compares well 
with the other numerical data as well as with the experimental results. For this example, it is 
also observed that the numerical simulation data do not match well with the experimental results 
for grazing incidence. 
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Figure 27 Thin, metallic diamond-shaped flat plate with a = 3.592 cm 
and descritized into quadrilateral patches 



Figure 28 Monostatic RCS of diamond-shaped plate (shown in Figure 27) when 

illuminated by H-polarized plane wave cTo = 0° at = 90° and = -90° , 

solid line- present method, short dashed line- RWG functions, dashed dot line 
ESP-code, solid square symbols- measurement, frequency = 1 1.81 IGHz. 
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Figure 29 Monostatic RCS of diamond-shaped plate (shown in Figure 27) when illuminated 
by E-polarized plane wave = 90° at = 90° and (j)- = -90“ , solid line- 
present method, short dashed line- RWG functions, dashed dot line ESP-code, 
solid square symbols- measurement, frequency = 1 1 .81 IGHz. 





Figure 30 Monostatic RCS of diamond-shaped plate (shown in Figure 27) when 
illuminated by H-polarized plane wave = 0° at = 0° and =180°, solid line- 
present method, short dashed line- RWG functions, dashed dot line ESP-code, 
frequency = 1 1 .8 1 IGHz. 
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Figure 31 Monostatic RCS of diamond-shaped plate (shown in Figure 27) when 
illuminated by E-polarized plane wave = 90"" at </>. = O'" and = 180^, solid 

line- present method, short dashed line- RWG functions, dashed dot line ESP- 
code, solid square symbols- measurement, frequency = 1 1.81 1 GHz. 
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Figure 32 Monostatic RCS of diamond-shaped plate (shown in Figure 27) when illuminated 
by E-polarized plane wave = 90" as a function of for frequency = 11.811 

GHz, 6 I = 90" . solid line- present method, short dashed line- RWG functions, 

dashed dot line ESP-code, solid square symbols- measurement, frequency = 
11.811GHz. 
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4 CONCLUSION 


The electric field integral equation approach in conjunction with the method of moments 
with quadrilateral meshing has been described to estimate the electromagnetic scattering from 
irregularly shaped thin metallic plates. Discretizing an irregularly shaped thin plate into 
quadrilateral patches and using proper expansion and testing functions over the quadrilateral sub- 
domains, the integral equation has been reduced to a matrix equation. From the knowledge of 
the surface currents density, obtained after solving the matrix equation, the EM scattering from 
various irregularly shaped plates has been determined and compared with earlier published 
results. The novelty in the present approach is the use of quadrilateral patches instead of well 
known and often used triangular patches. It has been observed that the use of quadrilateral 
patches instead of triangular patches reduces the number of unknowns and hence it requires 
handling of reduced size matrices. Also, the quadrilateral patches are more natural to satisfy the 
boundary edge currents condition, i.e., zero normal component compared to the triangular 
patches. This property makes it possible to achieve convergent results fewer quadrilateral 
patches. 

The numerical results obtained using quadrilateral patches have been compared with the 
numerical results obtained using the triangular patches and ESP-code. The numerical results 
obtained by the present code have also been compared with the measured data obtained in the 
NASA Langley Experimental Test Range Facility. 


APPENDIX A 

Construction of the vector basis function R„, (x,y) for an arbitrarily shaped quadrilateral 
in the x-y coordinate system is a complex problem. However, a coordinate transformation that 
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transforms a quadrilateral in the x-y plane into a square element in a new (^, rj) coordinates 
plane simplifies the derivation as shown below. Using the following transformation 


where 


^ and y = N. (^, 1^) y. 

/=1 /=! 


(Al) 


N, =l(l-^)(l-; 7 ), N, + , N, =l(l + ^)(l + ; 7 ) , and N, =^(l + a(l-7) 

4 4 4 4 

With the coordinate transformation given in (Al ), the n'’' quadrilateral is transformed into a 

square with sides equal to 2, as shown in Figure 3. The vector basis function (V, along the edges 

of the quadrilateral can be shown to be [ 6] 


yv, =-^(l-aV7, A, =^(1 + /7)V^, N,=-Ml + ^)Vri , (A2) 

4 4 ■ 4 4 

If h is the unit normal vector to the n'*' quadrilateral, then the basis function representing the 

surface current can be written as 




(A3) 



(A4) 


where the various quantities appearing in the above expressions can be expressed as [6 ] 


V^ = (x|^-5)|^)/det[y] , 
hri hri 


V;? = (-x|| + y||)/det[y] 

a^ 


det[7} = 


dx dy Bx By 


(A5) 


(A6) 


With simple mathematical manipulation it can be shown that 
4 det[7] 


(A7) 
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